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Jackiw-Nohl-Rebbi two-instanton as a source of magnetic monopole loop 
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We demonstrate that a Jackiw-Nohl-Rebbi solution as the most general two-instanton generates 
a circular loop ol magnetic monopole in four-dimensional Euclidean SU(2) Yang-Mills theory, in 
contrast to the one-instanton solution in the regular gauge for which no such magnetic monopole 
loops exist. These results together with our previous result indicate that two-instanton solution 
and two-meron solution with the same asymptotic behavior in the long distance are responsible for 
quark confinement based on the dual superconductivity picture. 

PACS numbers: 12.38.Aw, 21.65.Qr 



I. INTRODUCTION 

The quark confinement criterion a la Wilson [l|l is given 
by the area decay of the Wilson loop average (VFci-^^]), 
i.e., the vacuum expectation value of the Wilson loop op- 
erator Wc[s^] written in terms of the Yang-Mills field 
[3] £/fj,{x). In order to clarify the mechanism of quark 
confinement, it would be efficient to construct ensembles 
of gauge field which exhibit confinement. One antici- 
pates that a class of solutions of the classical Yang-Mills 
equation may give the dominant contribution to quark 
confinement. Among them, especially, topological (soli- 
ton) solutions [j] with non-trivial topological charge Qp, 
i.e., the Pontryagin index, are good candidate for a first 
examination. 

The well-known topological solutions are instantons 
pl-[lO| and merons |lll - [l3| : The instanton is a solution 
of the self-duality equation *F — zLF (first order dif- 
ferential equation) and has an integer-valued topologi- 
cal charge Qp = ±1, ±2, • • • , while the meron is a solu- 
tion of the original second order differential equation ob- 
tained without imposing the self-duality condition and 
has a half-integer- or integer-valued topological charge 
Qp = ±1/2, ±1,---. Remarkably, the instanton has a 
finite Euclidean action proportional to the topological 
charge S — 8Tr'^\Qp\/g'^ due to self-duality, while the 
meron has the logarithmically divergent Euclidean action 
due to short distance singularity, since the meron has no 
collective coordinate corresponding to the size. There- 
fore, the requirement of finite action inevitably excludes 
the meron solution from the candidates. However, it has 
been pointed out in [iJl that merons can be responsible 
for the occurrence of the confinement phase, if we take 
into account the free energy based on the action-entropy 
argument, since the logarithmic divergent action is com- 
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parable with the entropy associated with the meron con- 
figurations calculable from the integration measure in the 
path-integral formulation. 

It is believed that a promising mechanism for quark 
confinement is a dual superconductivity [15|, |la| . For this 
mechanism to work, however, magnetic monopoles must 
exist and be condensed to cause the dual superconductiv- 
ity. In this picture, magnetic monopoles are considered 
to be the most important degrees of freedom relevant to 
confinement in the dual description. Therefore, we are 
lead to estimate how magnetic monopoles contribute to 
the Wilson loop average to fulfill the confinement crite- 
rion [17|. 

On the other hands, we can ask which configuration 
of the Yang-Mills field can be the source for such mag- 
netic monopoles relevant to confinement in the dual de- 
scription. The simplest configuration examined first was 
the one instanton solution [j|. However, it has been 
confirmed in Jl^, |2C)| - [22| that magnetic monopole loops 
are not generated from one-instanton configuration, as 
briefly reviewed in the Introduction of our previous pa- 
per [l^|. Note that the magnetic monopole expressed by 
the current fc is a topological object of co-dimension 3, 
therefore, it is a one-dimensional object (a closed current 
due to the topological conservation Sk = 0) in the four 
dimensional space, while it is a point-like object in the 
three dimensional space. 

In previous papers, we have shown that two meron 
solution with a unit total topological charge \Qp\ = 1 
leads to circular loops of magnetic monopole joining a 
pair of merons in an analytical way [l8| and a numerical 
way [l9j, although one-instanton and two-merons have 
the same total topological charge. This result is in good 
agreement with the numerical result [23j obtained in the 
Laplacian Abelian gauge. 

In this paper we examine the two-instanton solution of 
Jackiw-Nohl-Rebbi (JNR) 9] from the viewpoint raised 
above. In the conventional studies on quark confine- 
ment, the multi-instanton solution of 't Hooft type [6| 
has been used extensively to see the interplay between in- 
stantons and magnetic monopoles |2lL l25l - l27| . However, 
the 't Hooft instanton is not the most general instan- 



ton solutions except for the one-instanton case in which 
the 't Hooft one-instanton agrees with the well-known 
one-instanton solution in the singular gauge. In con- 
trast, the JNR two-instanton solution is the most general 
two-instanton solution with the full collective coordinates 
(moduli parameters), while the 't Hooft two-instanton 
solution [611 is obtained as a special limit of the JNR so- 
lution. We demonstrate in a numerical way that a cir- 
cular loop of magnetic current k is generated for a JNR 
two-instanton solution. In addition, we present the con- 
figuration of the color field which plays the crucial role 
in our formulation. The implications of this result for 
quark confinement will be discussed in the final section. 
Incidentally, the JNR two instanton was used to study 
the relationship between dyonic instantons as a super- 
tube connecting two parallel D4-branes and the mag- 
netic monopole string loop as the supertube cross-section 
in (4+1) dimensional Yang-Mills-Higgs theory f28|, since 
dyonic instantons of 't Hooft type do not show magnetic 
string and D4-branes meet on isolated points, instead of 
some loop. These facts became one of the motivations to 
study the JNR solution from our point of view. 



II. MAGNETIC MONOPOLE CURRENT ON A 
LATTICE 

It is a nontrivial question how to realize the magnetic 
monopole in the pure Yang-Mills theory without matter 
field, while in the Yang-Mills-Higgs system such as the 
Georgi-Glashow model, the magnetic monopole in a non- 
Abelian gauge theory has been constructed long ago ^3] . 

In pure Yang-Mills theory in the absence of matter 
fields, two methods are currently known (to the best 
of our knowledge) for extracting magnetic monopole de- 
grees of freedom: 

1. Abelian projection due to 't Hooft [30] 

2. Field decomposition due to Cho, Duan and Ge, 

Hi-Ill 



Faddeev and Niemi, Shabanov 



The first method (Abelian projection) has succeeded 
to exhibit the Abelian dominance [S^, [SJ] and mag- 
netic monopole dominance [33| for quark confinement by 
adop ting the Maximally Alaelian (MA) gauge [3l|. See 
e.g., [35| for a review. However, the following questions 
were raised for results obtained in MA gauge, which are 
fundamental questions to be answered to establish the 
gauge-independent dual superconductivity in Yang-Mills 
theory. (1) How to extract the "Abelian" part respon- 
sible for quark confinement from the non- Abelian gauge 
theory in the gauge-independent way. (2) How to define 
the magnetic monopole to be condensed in Yang-Mills 
theory in the gauge-invariant way even in absence of any 
fundamental scalar field, in sharp contrast to the Gcorgi- 
Glashow model. 

The MA gauge is reproduced as a special limit of 
the second method. In fact, the first method is noth- 
ing but a gauge-fixed version of the second method. In 



other words, the second method is a manifestly gauge- 
covariant reformulation of the first one. Therefore, the 
second method enables us to answer the above questions. 
From this viewpoint, the second method has been devel- 
oped in a series of our papers J40h42| . Ilj - lijl . for SU{2) 
gauge group and [4J, |43, Ho] for SU{3) or SU{N) gauge 
group. In the second method, a magnetic monopole can 
be defined in a manifestly gauge-invariant way using new 
variables obtained from the original Yang-Mills field by 
change of variables. Moreover, the Wilson loop opera- 
tor can be exactly rewritten in terms of the magnetic 
monopole defined in this way through a non-Abelian 
Stokes theorem [5ll - l53j . 

In order to define the magnetic monopole on a lattice, 
we recall the second method based on a non-linear change 
of variables in a continuum SU{2) Yang- Mills theory. We 
introduce a color field n(a;) with a unit length: 

nix)=n^ix)T'', (T^ := a^/2) (1) 

n^ix)n'^ix)^l, (2) 

where aA {A — 1, 2, 3) are Pauli matrices. 

The color field is determined by imposing a condition 
which we call the reduction condition. A reduction con- 
dition is given by minimizing the functional 



Fred := f d^x^tT[{D^[A]n{x)}% 



(3) 



The local minima are given by the differential equation 
which we call the reduction differential equation (RDE) 



- D^[A]D^[A]n{x) ^ \{x)n{x). 
In this theory, a composite field 



(4) 



V^(a;) := Cf,{x)n{x) ~ ig ^ [df,n{x),n{x)] (5) 

play an important role where c^(a;) :— 2tr (n(a:;)A^(a;)). 
For instance, the field strength of V^(a;) is parallel to 
n(a;): 

F^.[V] = d^Y, - d,Y^ - ig [V^, V,] 

= {dt_,Ci, - dyC^ + 2ig'^iv{n [d^n, dun\ ) } n 
:= G^.n. (6) 

Thus, we can define the gauge-invariant field strength 
G^iy[x) = 2tr(nF^i,[V]) and the gauge invariant 
monopole current as 

A:^(x) := a/G^''(a;) = \e^''""' d,Gp„{x). (7) 

Once the reduction condition is solved, thus, we can ob- 
tain the monopole current fc^ from the original gauge 
field Ap. The gauge- invariant magnetic charge is defined 



g™ := / (fa^kf^ix), 



(8) 



in a Lorentz (or Euclidean rotation) invariant way |T7| . 

In this paper, we carry out the procedures explained in 
the above in a numerical way. We use the lattice regular- 
ization for numerical calculations where the link variable 
Ux,fi. is related to a gauge field in a continuum theory by 

{I'X+ajl \ 

ig / dyA^{y) I , (9) 

where P represents a path-ordered product, a is a lattice 
spacing and jl represents the unit vector in the ^ direc- 
tion. The lattice version of the reduction functional in 
SU{2) Yang-Mills theory is given by 

F,.od[n, [/] = ^ {1 - 4tr (C/,,^n,+,^C/]_^n,) /tr (1)} , 

(10) 



x,n 



where n^ is a unit color field on a site x, 

nx=n^T^, nini = l. (11) 

We introduce the Lagrange multiplier A^; to incorpo- 
rate the constraint of unit length for the color field ([TT]) . 
Then the stationary condition for the reduction func- 
tional is given by 



d 



dr 



|F„d[n,L/]-i^A,(n^n^- 1)1=0. (12) 



When Fred takes a local minimum for a given and fixed 
configurations {Ux^fi}, therefore, a Lagrange multiplier 
Xx satisfies 



and the color field n^ satisfies 
where 



(13) 



(14) 



W, 



-/=4^tr(L/.,^n.+aAf/i,^r^ 



M=l 



Ux-af..,T^Ul_^.^^nx-af.]/tr{l). (15) 



Eq. ([T3| is a lattice version of the reduction differential 
equation (RDE). We are able to eliminate the Lagrange 
multiplier to rewrite (1131) into 



W' 



' VW^W 



(16) 



A derivation of this equation is given in Appendix A. The 
color field configurations {n^;} are obtained by solving 
(1161) in a numerical way. 

After obtaining the {n^;} configuration for given con- 
figurations {Ux,fj.} in this way, we introduce a new link 



variable 14.^, on a lattice corresponding to the gauge po- 
tential (P iDy 



T/ ^'■^,fJ. 

yx,tl — r- J 

-'-^x.fi •' ^x,^. I ^x'^ x.ii^x-\-a(l- 



(17) 



Finally, the monopole current k^,^ on a lattice is con- 
structed as 



v. p. (7 



4tt 



through the angle variable of the plaquette variable 

^x,Mn,V] 
= a-2 arg (tr {(1 + n,) Vx,pVx+af.,uV^+ao,pVl,] 

/tr(l)). (19) 

In this definition, kx.p. takes an integer value |44l.l4a|. 

To obtain the {nx} configuration satisfying ((T5)) . we 
recursively apply (|16p to xix on each site x and update it 
keeping n^; fixed at a boundary dV of a finite lattice V 
until Frod converges. Since we calculate the {kx,p,} config- 
uration for the instanton configuration in this paper, we 
need to decide a boundary condition of the {nx} configu- 
ration in the instanton case. We recall that the instanton 
configuration approaches a pure gauge at infinity: 



gA^ix) ^ th\x)d^h{x) + Oi\x\-^). 



(20) 



Then, n(a;) as a solution of the reduction condition is 
supposed to behave asymptotically 

n{x)^h^x)T3h{x) + Oi\x\-''), (21) 

for a certain value of a > 0. Under this idea, we adopt a 
boundary condition as 

n^™"'^ := h\x)T3h{x), x € dV. (22) 

In practice, we start with an initial state of the {n^;} 
configuration: nj^f'* = h^ {x)T3h{x) for x £ V. Then, 
we repeat updating n^; on each site x according to (J16p 
except for the configuration n^°""'^ on the boundary dV. 
It should be remarked that these asymptotic forms ([20]) 
and (PT|) satisfy the RDE asymptotically in the sense that 



D^[A]nix) -^ (l^l -^ oo) 



together with 



A(a;)^0 (|x|^oo) 



(23) 



(24) 



which is necessary to obtain a finite value for the reduc- 
tion functional ]1§\ 



Fred = / d'''x-X{x) < CX). 



(25) 



p=5, b,=(0, 0, 0, 0) 



p=10, b,=(0, 0, 0, 0) 





(a) 



(b) 



p=15, b,=(0, 0, 0, 0) 



p=20, b„=(0, 0, 0, 0) 





(c) 



(d) 



FIG. 1: One instanton in the regular gauge and the associated magnetic-monopole current fc^.^ for various choice of size 
parameter p: (a) p = 5a, (b) p = 10a, (c) p = 15a and (d) p = 20a. The grid shows an instanton charge density D^ on X\-X2 
{xz = Xi = Q) plane. The black line on the base shows the magnetic monopole loop projected and colored lines shows a contour 
plot of the instanton charge density. Figures are drawn in units of a. 



III. ONE INSTANTON IN THE REGULAR 
GAUGE 



where 1 is a 2 x 2 unit matrix, and 



The one-instanton solution in the regular (or nonsin- 
gular) gauge is specified by a constant four-vector repre- 
senting the center (6^, b^,b'^,b*) e M^ and a positive real 
constant representing the size (width) p > 0: 






(26) 



h''{x)T3h{x) 



2(2;ia;3 - 2:2X4) 



Ti 



2(xi.X4 + X2X3) 



To 



'-T, 



(29) 



where |a:p = Xfj_Xp, is the standard Euclidean norm and 
7/|Ui) ' is the symbol defined by 



lu!)'^^ = ^Afi^4 ± SAfj.5^4 T Sa^Si 



p4- 



In this case, we obtain from pO|) and ((2T|) : 



h{x) = 7-76^, (Cfj, = {-iat, 1)) , 



(27) 



(28) 



This exactly agrees with the standard Hopf map [5J| . 

The topological charge density is maximal at the point 
X — b and decreases algebraically with the distance from 
this point in such a way that the instanton charge Qv 
inside the finite lattice V = [—aL, aL]'^ reproduces the 
total instanton charge Qp = 1. We construct the instan- 
ton charge Qv on a lattice from the configuration of link 



variables {Ux.fj.} according to 

Qv =a'^ ^ D^, 
xe{v-dv} 

,-4 ±4 



(30) 



D,::= 



1 a- 



E 



hVyi. '~^x,^v^x,p<T }■> 



fi.i^,p,(T—:Ll 



u. 



'^ x.ii^ x-\-aLi,i/^ r 



x+aO-ii 



^ x.v 



(31) 
(32) 



where Dx is a lattice version of the instanton charge den- 
sity, and V — dV represents the volume without a bound- 
ary, and e is related to the usual e tensor by 



gM!^P<T .^ sgn(/x)sgn(z^)sgn(p)sgn(cr)el^ll'^ll''ll'"l. 



sgn(Ai) := -TT- 

ImI 



(33) 



Our interest is the support of fc^./i, namely, a set of 
links {x, /i} on which fc^j^^ takes non-zero values k^^^ ^ 0. 
This expresses the location of the magnetic monopole 
current generated for a given instanton configuration. By 
definition (J18p . the number of configurations {fc^.^j} are 
(2L)'* X 4. (The number of configurations {kx,^i} is not 
equal to (2L + 1)* x 4, because we can not calculate kx^^j. 
at positive sides of the boundary dV due to the definition 
of (fTS]) based on the forward lattice derivative.) 
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\kx,p,\ > 
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fco; 


M — 


-1 k.,^ = Q 


k. 


,M = 1 


Qv 
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59105336 
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0.9674 


10 









8 


59105328 
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0.9804 


15 









12 


59105320 




12 


0.9490 


20 









12 


59105320 




12 


0.8836 



TABLE I: The distribution of kx,^i, and the instanton charge 
Qv for a given r. 

The results are summarized in TABLE|T]and FIG.[T] In 
our calculations of the monopole current configuration, 
we fix the center on the origin 



(6i,62,&3^&4) = (0,0,0,0), 



(34) 



and change the value of p. 

For a choice of L = 31, the total number of config- 
urations {kx,^,} are 62* x 4 = 59105344. Although the 
current kx,ii is zero on almost all the links (a:,/i), it has 



a non-zero value \kx^pi 
e.g., 4-1-4 links for p 



= 1 on a small number of links, 
5a. The number of links with 



-1 is equal to one with kx.. 



-1, which re- 



flects the fact that the current kx.,_i draws a closed path 
of links. It turns out that there are no configurations 
such that \kx,^\ > 1. 

We see how Qv reproduces the total instanton charge 
Qp — 1 for the choice of p. The instanton charge den- 
sity is equal to zero when A^(x) is a pure gauge, i.e., 
F^i/(x) = 0. Therefore, the more rapidly A^(a;) con- 
verge to a pure gauge, the more precisely Qy reproduces 



the proper value Qp for a given L. For one-instanton in 
the regular gauge, clearly, A^(a;) more rapidly converge 
to a pure gauge for smaller p. Indeed, Qv for p = 10a 
reproduces Qp = 1 more precisely than p = 15a and 
p — 20a. On the contrary, Qy for p = 5a reproduces 
Qp = 1 poorly compared with p — 10a. This is be- 
cause the lattice is too coarse (a lattice spacing is not 
sufficiently small) to estimate properly the rapid change 
of the instanton charge distribution concentrated to the 
neighborhood of the origin for small p. 

If A^(a;) converge rapidly to a pure gauge ih){x)d^h{x) 
so that Qv gives a good approximation for an expected 
value Qp, then the color field is expected to behave as 
n(a;) — >■ K^{x)T^h(x) asymptotically and our choice of 
the boundary condition n^°""'^ := h\x)T^h{x) at the 
boundary x G dV is well motivated. 

In FIG. [1] the support of fc^,^ is drawn by projecting 
the four-dimensional space on the 14 = hyperplane 
(3-dimensional space) for the choice of p = 5a, 10a, 15a 
and 20a. This figure shows that the non-zero monopole 
current forms a small loop. The size of the magnetic 
monopole loop hardly change while p increase. This is 
an indication that the magnetic monopole loop for one- 
instanton solution disappears in the continuum limit of 
the lattice spacing a going to zero. 



IV. TWO INSTANTON OF THE 
JACKIW-NOHL-REBBI (JNR) TYPE 

It is known by the ADHM construction [l3| that the N- 
instanton moduli space has dimension 8iV. For TV = 1, 
the 8 moduli parameters are interpreted as 4-f 1-1-3 de- 
grees of freedom for the position, size and SU(2) orienta- 
tion (global gauge rotations), respectively. 

For A^ = 2, the Jackiw-Nohl-Rebbi (JNR) instanton [i| 
is the most general charge 2 instanton as explained below. 
The explicit form of the JNR two-instanton solution is 



5A^(x)=-TV^^a.ln0, 



JNR 



'-^ '1,1V VmK / . 



2pl [x'' - 6^) 

^0 i\^-br?? ' 



OJNR • = 



E 



pI 



6.P' 



(35) 
(36) 

(37) 



The JNR two-instanton has 4x3 + 3 + 3 = 18 
parameters, which consist of three pole positions 
{bl bl bl bt), [blbl blbf), (blblblbj) and three scale 
parameters PQ,pi,p2 including the overall SU{2) orien- 
tation. Note that the number of poles is one greater 
than the number of the instanton charge. Although the 
parameter count of the JNR two-instanton appears to 
exceed the 16 dimensions of the N = 2 moduli space, the 
JNR two-instanton has precisely the required number of 
parameters for the N = 2 general solution. In fact, one 
parameter is reduced by noting that the multiplication 
of the scale parameter by a constant does not alter the 
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FIG. 2: JNR two-instanton and the associated circular loop of the magnetic monopole current kx,^- The JNR two-instanton 
is defined by fixing three scales po ~ pi ~ P2 = 3a and three pole positions &q,&j,62 which are arranged to be three vertices 
of an equilateral triangle specified by r: (a) r = 5a, (b) r — 10a, (c) r — 15a and (d) r = 20a. The grid shows an instanton 
charge density D^ on a;i-a;2 (a;3 = a;4 = 0) plane. The associated circular loop of the magnetic monopole current is located on 
the same plane as that specified by three poles. The black line on the base shows the magnetic monopole loop projected on the 
2:1-2:2 plane and the arrow indicates the direction of the monopole current, while colored lines on the base show the contour 
plot for the equi-D^ lines. Figures are drawn in units of a. 



solution, so only the ratios Pr/po ('^ = lj2) are relevant. 
Moreover, one of the degrees of freedom corresponds to 
a gauge transformation Q . 

The 't Hooft two-instanton which is more popular and 
has been used extensively in the preceding investigations 
is given by 



OtHoott 



2 



(38) 



The 't Hooft two-instanton has only 4 x 2 -f 2 -f 
3 = 13 parameters, which consist of two pole posi- 
tions (6j, 6^, 6f , &f), (62J &2: ^2: ^1) ^^^^ ^wo scalc param- 
eters pi,p2 including the overall SU{2) orientation. The 
't Hooft solution is reproduced from the JNR solution in 
the limit po = \bQ\ ^ 00, namely, the location of the first 
pole 60 is sent to infinity keeping the relation po = l^ol- 

The crucial difference between the JNR and the 't 
Hooft solutions is the asymptotic behavior. The JNR 
solution goes to zero slowly, while the gauge potential 



produced by the 't Hooft ansatz tends rapidly to zero at 
spatial infinity |a;| — ?► cx). In fact, for iV = 1 the 't Hooft 
ansatz gives the one-instanton in the singular gauge with 
the asymptotic behavior 



A^(x)- 0(1x1-3) \^\^^^ 



(39) 



while the one-instanton in the regular gauge exhibits the 
asymptotic behavior 



A^{x) ^ 0{\x\-^) |x|^oo. (40) 



In the case of the JNR two-instanton, we obtain 



Hx) = -rn-e/x' (e/x = '^l = (^"^^^ 1)) (41) 
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(a) 



(b) 



FIG. 3: The configuration of tlie color field w^ = (r 



and a circular loop of the magnetic monopole current kx 



obtained from the JNR two-instanton solution (r = 20a), viewed in (a) the x^-xz {x\ = a;4 = 0) plane which is off three poles, 
and (b) the x\-xz (x2 = X4 = 0) plane which goes through a pole bp. The magnetic monopole current k^^^j, and the three poles 
of the JNR solution are projected on the same plane. Here the SU{2) color field {nl.,n'?^, n^) is identified with a unit vector in 
the three-dimensional space (xi, 0:2,3:3). Figures are drawn in units of a. 



and 



htf \rr Ul \ "2 {X1X3 + X2X4:) 

h\x)T-ih{x) = 2 ^1 



2 {-XiXA + 0:2X3) 



x^ 
^2 



+ 



T2 

2 

J. 3 



(42) 



This is another form of the standard Hopf map [5J] . 

We equate three size parameters and put three pole po- 
sitions {blblblbt), Kblblbi), (blblblb^^) on the 
2:3 = 0:4 = plane, so that the three poles are located at 
the vertices of an equilateral triangle: 



po =pi = P2 = p, 
iblbl,blbt)=ir,0,0,0) + \ 



A, 



{blblblbt) 



-i-f-'°'» 



(43) 
(44) 

(45) 
(46) 



where A is a small parameter introduced to avoid the 
pole singularities aX x — br- 

In the numerical calculation, we choose 



(0.1a, 0.1a, 0.1a, 0.1a). 



(47) 



Then we have searched for monopole currents by chang- 
ing r. We have checked that for a given size of the lattice 



r 


\kx,^i\ > 1 k 


x,n — 


J- f^x,fi — U rCa 


,M = 1 


l/r 


Qv 


5 





14 


59105316 


14 


5.6 


1.903 


10 





26 


59105292 


26 


5.2 


1.969 


15 





40 


59105264 


40 


5.3 


1.950 


20 





51 


59105242 


51 


5.1 


1.862 



TABLE II: The distribution of kx,^ and the instanton charge 
Qv for a given r. 



L there is a suitable range of p, in which the effect of 
the discretization is sufficiently small. In view of this, we 
adopt L = 31 and p = 3a. 

The results are summarized in TABLE HIl and FIG. [H 
and |3l In FIG. [U as in the case of one-instanton in the 
regular gauge, we draw the distribution of the instan- 
ton charge density Dx and the support of the magnetic 
monopole current k^^^, projected on the 2:4 = hyper- 
space for r — 5a, 10a, 15a, 20a and p — 3a. 

The instanton charge density Dx of the JNR two- 
instanton takes the maximal value at a circle with radius 
i?/, rather than on the origin, on the a;i-a:2 plane. This is 
not the case for the 't Hooft two instanton in which the 
instanton charge distribution concentrates near the two 
pole positions, as is well known. 

We have found that non-vanishing monopole currents 
originating from the JNR two-instanton forms a circular 
loop. The circular loops of the magnetic monopole cur- 
rent are located on the same plane as that specified by 
three poles 60,61,62- The size of the circular loop, e.g., 



8 



the radius R, increases proportionally as r increases and 
the circular loops constitute concentric circles with the 
center at the origin, within the accuracy of our numerical 
calculations. 

To reproduce the correct distribution of the instanton 
charge density Dx in the numerical calculations, we need 
to further approximate the link variable Ux,^, defined by 
^ using a discretization with a width e = a/n for the 
line integral on a link: 



Ux., 



exp < le 



n 



k=l 



X + {k — l)e/i) 



-gAfj^ix + kefj.)] 



(48) 



We have chosen the lattice spacing a — 1 and the total 
number of partition points n — 20. If n is relatively small, 
the numerical result for the instanton charge density D^ 
gives a wrong distribution with remnant of three poles. 

In TABLEini we find that the ratio between the length 
£ := ^^ \kx,p,\ of the magnetic monopole current kx.f_i 
and the size r of the equilateral triangle of JNR is nearly 
constant. 



which implies 



/r~5.2, 



R^/r ~ 0.65, 



(49) 



(50) 



where we have used a relation £ ~ 8Rm for large £ since 
a closed current kx,ij, consists of links on a lattice and 
the relation £ ~ 2'iTRm in the continuum does not hold. 
Moreover, the magnetic monopole loop passes along the 
neighborhood of contour giving the absolute maxima of 
the instanton charge density. 



i?„/i?/ ~ 0.65/0.54 ~ 1.2, 



(51) 



where Rj ~ 0.54r. 

FIG. [3] shows the relationship between the magnetic 
monopole loop kx^^ and the color field ria; configuration. 
The vector field {na;} is winding around the loop, and it is 
indeterminate at points where the loop pass. The config- 
urations of the color field giving the magnetic monopole 
loop were made available for the first time in this study 
based on the new reformulation of Yang-Mills theory. 



we have constructed the magnetic monopole current k^ 
on a dual lattice where the resulting magnetic charge is 
gauge invariant and quantized according to the quanti- 
zation condition of the Dirac type. 

For the two-instanton solution of the Jackiw-Nohl- 
Rebbi type, we have discovered that the magnetic 
monopole current A:^ has the support on a circular loop 
which is located near the maxima of the instanton charge 
density. Thus, we have shown that the two-instanton so- 
lution of the Jackiw-Nohl-Rebbi type generates the mag- 
netic monopole loop in four-dimensional SU(2) Yang- 
Mills theory. In the same setting, we have found that 
the magnetic current has the support only on a plaque- 
tte around the center of the one-instanton. This result 
confirms that no magnetic monopole loop is generated 
for one-instanton solution in the continuum limit. 

Combining the result in this paper with the previous 
one [18] , we have found that both the JNR two-instanton 
solution and two-mcrons solution with the same asymp- 
totic behavior at spacial infinity: 



A^{x) ^ 0{\x\-^) |a;|^oo 



(52) 



generate circular loops of magnetic monopole, which 
should be compared with the 't Hooft (multi) instanton 
with the asymptotic behavior at spacial infinity: 



A^{x) ^ 0{\x\'^) \x\ 



— > oo. 



(53) 



We expect that these loops of magnetic monopole are 
responsible for confinement in the dual superconductivity 
picture. This result seems to be consistent with the claim 
made in [53,|5y|. However, the correspondence between 
the instanton and magnetic monopole loop is not one-to- 
one. To draw the final conclusion, we need to collect more 
data for supporting this claim. Moreover, it is not yet 
clear which relationship between instanton charge and 
the magnetic charge holds in our case, as studied in [53- 
'59*1 . These issues will be investigated in future works. 
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V. CONCLUSION AND DISCUSSION 

For given instanton solutions of the classical Yang- 
Mills equation in the four-dimensional Euclidean space, 
we have solved in a numerical way the reduction condi- 
tion to obtain the color field which plays the key role 
to define a gauge-invariant magnetic monopole in our re- 
formulation of the Yang-Mills theory written in terms 
of new variables. Here we have used a lattice regular- 
ization |44l | for performing numerical calculations. Then 



Appendix A: Deriving the reduction equation on a 
lattice 

Taking the square root of both sides of P^ yields 

A, = ±^Jw^W^, (Al) 

which is substituted into ([T3l) to obtain 



n^ = ±- 



^jWwi 



(A2) 



We must choose a correct sign in the right-hand side of correct sign of the right-hand side by examining whether 

the equation (|A2|) so that a solution of this equation gives the equation (|A2p converges to the RDE in the contin- 

a solution of the original equation ([T3l) . uum limit a -^ 0. By substituting Ux^fj. = e"'"^^'''^^) into 

In what follows, we adopt a direct method of choosing a (jlSp and expanding it in powers of a, W,^ is written as 






4^tr 



a 



n{x + ap.) - ia [gA^(a;), n(x -f- a/i)] — [gAf^{x), [gA^{x), n{x + a/i)] 

+ n{x — afi) + ia [gAfj_{x — afi), n{x — afi)] 



[gA^{x - a/i), [gAf,{x - a/i), n{x - a/i)]] ^ T 



/tr(l) + 0(a3) 



= 4 ^ tr [{2n(a;) + 0^9^9^11(0;) 

- 2ia^ [gA^(a;), 9^n(a;)] - ia^ [5apAp(x), n(a;)] 

-a2 [gA^(a;), [gA^(x), n(x-)]]} T^] /tr (1) -t- ©(a^) 
4 
= 4^ tr [{2n(x) + a2i^^[A]i?^[A]n(a;)} T^] /tr (1) + 0(a3) 

4 

= 8n^(x) + a' ^ {D^[A]D^[A]n{x)f + 0{a^). 
Then, the magnitude of W^ is given by 



'W^W^ = /64 + 16a2 ^ n'4(x) (L>^[A]i:'^[A]n(.T))-^ + 0{a^) 
Substituting (jA3l) and (|A4[) into (jA2p . therefore, we obtain 



(A3) 



(A4) 



8(-n^(x)±n^(x))+aM -^(i?,,[A]I?^[A]n(a;))^±n^(a;)^n^(x)(i^^[A]i?,JA]n(x))^ +O(a3) = 0, (A5) 



where the double-signs ± in (jA2p and (jA5l) correspond 
to each other. If we choose the minus sign, we have an 
unrealistic result n (x) = 0. Thus, by choosing the plus 
sign in (|16p . the leading term vanishes and the next-to- 
leading term leads to 

4 
J2 {D,[A]D,[A]n{x)f = -X{x)n^{x) (A6) 

tj.=i 

X{x) = - ^ n^{x) {D^[A]D^[A]n{x)f , (A7) 



r 



which is nothing but the RDE in the continuum theory: 



D^[A]D,,[A]n{x) = X{x)n{x). 



(A8) 
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